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Motivated by recent results of A'^th order muffin-tin orbital (NMTO) implementation of the density functional 
theory (DFT), we re-examine low-temperature ground-state properties of the anti-ferromagnetic insulating phase 
of vanadium sesquioxide V2O3. The hopping matrix elements obtained by the NMTO-downfolding procedure 
differ significantly from those previously obtained in electronic structure calculations and imply that the in- 
plane hopping integrals are as important as the out-of-plane ones. We use the NMTO hopping matrix elements 
as input and perform a variational study of the ground state. We show that the formation of stable molecules 
throughout the crystal is not favorable in this case and that the experimentally observed magnetic structure 
can still be obtained in the atomic variational regime. However the resulting ground state (two t2g electrons 
occupying the degenerate Cg doublet) is in contrast with many well established experimental observations. We 
discuss the implications of this finding in the light of the non-local electronic correlations certainly present in 
this compound. 



PACS numbers: 71.27.H-a, 71.30.H-h, 75.25.H-z 

I. INTRODUCTION 

Vanadium sesquioxide V2O3 holds a a very special place 
among the variety of physical systems exhibiting a metal- 
insulator transition (MIT). This compound is considered as a 
prototype of a Mott-Hubbard system: it displays a MIT from 
a paramagnetic metallic (PM) phase to an antiferromagnetic 
insulating (AFI) phase at low temperatures (ss IbQK), and a 
transition from a PM phase to a paramagnetic insulating (PI) 
phase at a higher temperature (w 500i4r)ji2iM In order to un- 
derstand several of its properties, the knowledge of its crystal 
structure plays a fundamental role : in the paramagnetic phase 
V2O3 can be characterized by a corundum cell in which V 
ions are arranged in V-V pairs along the c-hexagonal axis and 
form a honeycomb lattice in the basal ah plane. All vanadium 
ions are equivalent: each V'^+— ion has 3fi^ configuration and 
is surrounded by a nearly perfect oxygen octahedron. How- 
ever, a slight displacement of vanadium ions away from the 
centers of their octahedra is at the origin of a small trigonal 
distortion which removes the three-fold degeneracy of the t2g 
manifold. A non-degenerate aig orbital and doubly degener- 
ate eg orbitals are then separated by an energy gap measur- 
ing the trigonal distortion. At. A first order structural transi- 
tion takes place upon cooling, and the system becomes mono- 
clinic. Simultaneously, a peculiar antiferromagnetic (AF) spin 
order emerges with ferromagnetically ordered planes antifer- 
romagnetically stacked perpendicularly to the monoclinic hm 
axis. 

In late seventies, Castellani, Natoli and Ranninger^ (CNR) 
developed the first realistic description of V2O3. They re- 
alized that the peculiar structure observed in the AFI phase 
could not be explained in terms of a single band Hubbard 
model, and that the introduction of the orbital degrees of 
freedom into the model was a necessary ingredient in order 



to explain the experimental findings. The CNR model has 
been considered as a reliable model of V2O3 until the end 
of the nineties, when several experiments, e.g., x-ray absorp- 
tion spectroscopy by Park et al,^ and non-resonant magnetic 
x-ray scattering by Paolasini et alj^^ have demonstrated that 
this model was in need of corrections. The failing of CNR 
model came from the underestimation of the value of Hund's 
coupling (J ^ 0.2 ~ 0.3 eV) and the subsequent result of 
a spin-state S" = 1/2 for 1/'^+— ions. Instead, all previous 
experiments independently showed that each vanadium ion 
has a spin S — 1, observation that called for a deep revi- 
sion of the theoretical description of V2O3. For this reason, 
many theoretical works followed i^''°'''i'^''^''^-'^ In particular, 
very recently, the nature of the MIT in V2O3 has been stud- 
ied by means of LDAh-DMFT approaches. Laad et al^^ have 
proposed that the MIT is of orbital-selective type, and that 
large changes in the orbital occupation at the transition oc- 
cur in response to the crystal field changes. Poteryaev et al^\ 
with a different LDAh-DMFT approach showed that the crys- 
tal field, the correlations, and the orbital degrees of freedom 
are strongly interrelated. They also proposed that the trigo- 
nal distortion, strongly enhanced by correlations, acts as an 
external field in the orbital Hilbert space, in the same way as 
a magnetic field acts in the spin space. These two work a'^''^ 
have suggested that drastic modifications in the hybridization 
of the orbitals due to the correlation effects and the trigonal 
splitting were at the basis of the change in the population of 
the orbitals across the transition. 

Despite all this progress in the understanding of the MIT 
in V2O3, the controversy about the nature of the low- 
temperature antiferromagnetic phase remains. Specifically, 
there are contradictory results on how the two electrons are 
distributed within the t2g manifold, and whether or not orbital 
ordering takes place. The LDA+U calculation of Ezhov et 
al. suggested that the S = I state with no orbital degener- 
acy and a purely Cg occupancy is a possible candidate for the 
low-temperature antiferromagnetic phase. This idea, mutatis 
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mutandis, is present also in the LDA+DMFT calculations of 
Ref. I1I21I . Mila et alX on the contrary used an old idea by 
Allen— that magnetic and optical properties of all the phases 
of V2O3 show a loss of t^'^+-ion identity, and suggested that a 
good candidate for a ground state is a state in which vanadium 
ions have S — 1 and form vertical bonds (molecules) with the 
total spin Stot = 2. In our previous study", we showed us- 
ing a different correlated model that the molecular state can 
be stabilized throughout the crystal only if the correlation en- 
ergy within the molecule is big compared to the in-plane in- 
teraction energy. This gave a certain region in the parameter 
space where the molecular solution could be realized. Later, 
Tanaka'^ proposed that, in order to explain a large contribu- 
tion from the orbital magnetic moment, < L > ~— 0.5/i_B 
(the minus sign reminds that this contribution is opposite to 
the spin contribution 2 < S > ^ 1.7^b), observed in the 
non-resonant magnetic scattering experiments by Paolasini et 
alJ^, the orbital wave function had to be complex. He argued 
that spin-orbit coupling plays an important role in determining 
both the ground state and the low energy excitation spectrum 
ofVaOa. 

This controversy of theoretical results calls for a better 
quantitative description of the electronic structure of V2O3, 
which is a starting point for any model-based analysis. Re- 
cently, the NMTO method has been applied to study the elec- 
tronic structure of V203.'^ Surprisingly, after early works by 
CNR^ and by Mattheis,'"* this is the first serious attempt to 
model complicated electronic structure of V2O3 by an effec- 
tive single particle Hamiltonian. The NMTO method pro- 
vides a way to derive an effective Hamiltonian with only few 
energy-selective Wannier-like orbitals by integrating out all 
others, less important degrees of freedom. The result is an 
accurate tight-binding (TB) modeling of V2O3. The NMTO 
approach shows that oxygen contribution to hopping matrix 
elements and renormalization effects due to hopping paths via 
the Cg-tails (in addition to direct V-V t2g-hopping paths) are 
important and were only partially taken into account in Refs. 
lIsllTill . As one can see from Table I, the TB hopping param- 
eters obtained by the NMTO-downfolding technique differ 
significantly from the estimates based on the TB fitting per- 
formed on the LAPW band calculation of Mattheiss'** or by 
CNR. The NMTO analysis shows that the hopping through 
Eg-tails mainly influences the parameters in the basal plane. 
This is a bonding effect, which leads to an increase of the 
in-plane matrix elements. The other contributions via oxygen 
tails affect mainly the hopping along vertical bonds. This is an 
anti -bonding effect which significantly reduces aig — aig hop- 
ping matrix elements. The latter leads to the decrease of the 
correlation energy of vertical pairs, in this way reducing the 
stability of the molecular state. Simultaneous increase of the 
in-plane interaction energy favors uncorrected atomic sites, 
and in the end potentially breaks the molecular state. This 
new qualitative picture of the ground state of V2O3 in the pa- 
rameter space, suggested by NMTO-downfolding method,'^ 
calls for re-examination of the variational analysis of the ef- 
fective Hamiltonian in the antiferromagnetic insulating state, 
which is the main purpose of this paper 

The structure of the paper is as follows. We first present a 



brief outline of the results of our previous work, Ref. ifTll] and 
set the notations. We then minimize the effective Hamiltonian 
we derived using the values of the hopping parameters from 
Ref. [ 1 3] and performing a variational calculation with two 
kinds of variational functions: molecular and atomic. We dis- 
cuss the result in the light of the experimental facts and try to 
draw some conclusions for the physics of V2O3. 

Table I. Transfer integrals (eV) from TB calculations used 
by CNR^, from LAPW-calculations by Mattheiss'"^, and from 
NMTO-downfolding technique by Saha-Dasgupta^^. Nota- 
tions are explained in Table II in the Ref. iflHl . 
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0.2 
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0.05 


0.05 
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-0.23 


-0.27 


-0.26 



II. MODEL, FORMALISM AND NOTATIONS 

We assume that the low-temperature AFI phase of V2O3 
can be described by a super-exchange spin-orbital Hamilto- 
nian, iJsE, that can be derived perturbatively from a three-fold 
degenerate t2g Hubbard Hamiltonian, Hh, as done in Ref. 
1 1 1], to which we refer for further details. We can write the 
Hubbard Hamiltonian as: 

Hh = H, + , (1) 

where the kinetic term = Y^rj' E™™'^ *5;j"'c+„,^Cj'„,,^ 
includes a summation over nearest neighbor (nn) sites, over 
orbital (m,m' = 1,2,3) and spin (a =tii) indexes. Hu 
describes the on-site Coulomb interactions Ui (for electrons 
in the same orbital) and U2 (for electrons in different or- 
bitals) and Hund's coupling J. The hopping integrals t""," 
(to, to' — 1,2,3) can be expressed via a reduced set of pa- 
rameters: /i and p for out-of-plane hopping (ie, within the 
molecule), a, (3, a and r for in-plane hopping. In Table 
I we present the numerical values of hopping matrix ele- 
ments obtained by CNR^ through unrestricted Hartree-Fock 
calculations, by Mattheiss'"* throug a fitting procedure of the 
LAPW-calculations, and by Saha-Dasgupta et al'^ through the 
NMTO-downfolding technique. We shall later use these pa- 
rameters in the variational analysis where we compare various 
ground-states energies. For the Coulomb repulsion U2 and 
Hund's coupling J Ezhov et al}^^ and Mila et al."^ suggested 
J l.OeV ,1/2 ^ 2.5 eV. Recent optical studies of Qazilbash 
et a9^, however, yielded a smaller value for the Hund's cou- 
pling J = 0.5 eV. We have described already" how the range 
of these two parameters is subjected to large fluctuations in 
the literature. In what follows we shall fix C/2 ^ 2.5 eV and 
consider a range of J ~ 0.4 1.0 eV . 

We further assume, as found experimentally, that at each 
site two t2g-electrons are bound into a 5 = 1 state. We 
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therefore derive from iJn an effective Hamiltonian, Hcs = 
i?SE + ^^trig to describe the insulating phase of V2O3, where 
HsE and i?trig are given by 
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Here 5j = 1 is the spin at site j and rijma describes the oc- 
cupation of the m orbital on site j by an electron with spin 



(k) 

a, whereas the Olj are orbital exchange operators presented 



in Appendix C of Ref. ifTTIl . Multiplied by the correspond- 



ing prefactors 



or — - 



they de- 



4U2+iJ 12 [72+2,7 ' 

fine the effective value of the exchange that depends on the 
orbital occupation of the two sites i and j along the bond di- 
rection < i j >. The spin terms can be simply considered 
as projectors on the FM and AFM state. The term i?trig de- 
scribes the trigonal distortion that splits the three degener- 
ate two-electron states: legC^) = |0), joige^ = | — 1), 
jaigCp = |1). The energy splittings due to this distortion 
are defined as Ai = A2 = and A3 = At > for orbitals 
TO = 1,2 and 3, respectively. At is comparable in magni- 
tude to the hopping integrals, which allows us to treat the term 
Htiie on the same level as Hse- 



III. RESULTS 

In this section, we compare orbital and magnetic ground- 
state configurations of the effective Hamiltonian Hcs as ob- 
tained by a variational analysis with TB parameters from 
Mattheiss'"* and NMTO'^- models. 

The trial wave function can be most generally written as 



(4) 



where \4>n) describes the orbital part and ) refers to the 
spin part of the wave function on a site n. 

In the following, we shall use as a variational wave func- 
tion I'I'n) either an atomic state or a molecular state, with n 
labeling an atomic or a molecular site, respectively. In both 
cases, the average value of the Hcs over the corresponding 
state takes the form: 



(«'„|(*™|£reff|^m)|*n) = 



(5) 



This factorization is possible only in the mean field approxi- 
mation, in which single-site wave-functions, l^'n), otherwise 
entangled, factor into orbital and spin parts. By this 
we explicitly neglect any coupled spin-orbit fluctuations. The 
spin averaging is straightforward: for a ferromagnetic bond 
(5„ ■ Srn + 2)HF = 3 and (S„ ■ S^-1)hf = 0, while 



for an antiferromagnetic bond, (S'„ • Sm 
(5„ • S,n — 1)/7F = —2. Orbital averaging in Eq.(5]l, however, 
requires some algebra and is discussed in Ref. ifllll for both 
atomic and molecular states. In the following, we will first av- 
erage over all possible ordered magnetic structures, and then 
for each type of magnetic ordering find an orbital configura- 
tion that minimizes the total energy be means of a variational 
procedure. Note, that we limit the number of magnetically 
ordered structures only to those that can be realized on the 
corundum unit cell. 

As a preliminary step, before the full minimization of the 
variational energy (|5]l, we can compare the correlation energy 
of the ferromagnetic state of the vertical pair, 2p^/{U2 — J), 
with the superexchange energy in the basal plane, in order 
to have an idea about the relative order of magnitude of 
atomic and molecular energies. We remind that the corre- 
lation energy is defined as the difference between the ex- 
act ground state energy, Ey — —{p ~ p)"^ /{U2 ~ J), and 
the ground state energy in the Hartree-Fock approximation, 
Ey^ ~ [p^ + ^^)/{U2 — J)- The superexchange energy in the 
basal plane was approximated by (a^ + t^)/(C/2 — J) in Ref. 
ifTTll . where we assumed that the contributions from [3 and a 
hopping matrix elements were negligible, and should rather 
be written as {(3'^ + T^)/(t/2 - J) with the NMTO values of 
Table I, that shows that (3 is no more negligible, whereas a it 
is. 

There are two qualitatively different regimes of solutions. 
If the correlation energy is larger, the most appropriate varia- 
tional wave function for the whole H^s must be constructed 
in terms of molecular units with S^' = 2 and orbital pait \ip'^) 
given by 



\r±)ab=^{\±l)a\Q)b+\±l)bna) 



(6) 



where a and h define two sites of the vertical molecule, and 
±1 and denote the value of z component of the pseudospin 
operator, defining the orbital state. 

If, instead, the values of the exchange energy in the basal 
plane are larger than the correlation energy, the appropriate 
variational wave function is atomic-like and should be written 
as 

IC) = cos0,|O), + sin0,(cosi/;,|l), + sini/;,| - 1),) . (7) 

Contrary to the molecular case, this wave function allows all 
three states |0)i, and | — 1)^ to be present without any a 
priori restriction on their relative weight. The relative weight 
of these three states is then determined through the minimiza- 
tion procedure with respect to the variational parameters 9i 
and V'i- 

As quantitatively shown in Ref. ifTlll . Eqs. (6.19) and 
(6.20), in order to determine the nature (atomiclike, molec- 
ularlike) of the ground state, we need to compare the ab- 
solute value of the ratio of the molecular correlation energy 
(per atom) and the in-plane exchange energy, computed with 
Mattheiss parameters and new NMTO set of parameters. We 
obtain p/i/(Q;2 + r2) = 1.7 and 0.26 for Mattheiss and NMTO 
set, respectively. The strong reduction of this ratio for NMTO 
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set suggests that NMTO parameters do not favor the formation 
of the molecular state. 

There are few possible reasons for such large difference in 
the results of the two approaches. The discrepancy may be 
the consequence of the fact that NMTO method more accu- 
rately treats the effect of trigonal distortion on hopping ma- 
trix elements. In the approach of Mattheiss the influence of 
the trigonal distortion on matrix elements was assumed to 
be small and was neglected. As one can see from Fig. 14 of 
Ref. [131], where structures with varying amount of distortion 
have been generated, the aig — aig hopping along vertical 
bond, p, decrease significantly with the increase of the trig- 
onal distortion: NMTO estimates are p{At = 0) = —1.2 
eV, while p{At = 0.27eV) = -0.51 eV. The same figure 
gives much smaller estimate of variation of p for CNR model: 
p{At = 0) = -0.8 eV and p(A( = 0.27eV) = -0.7 eV. 
Other reasons of discrepancy, as pointed by the authors of 
Ref. 1 13] themselves, can be traced down to various approxi- 
mations used in constructing effective orbitals in the NMTO- 
downfolding, and to imprecise knowledge of the value of the 
covalency V(d)-0(p) mixing parameter and the charge trans- 
fer 3c? — 2p gap, extracted from nuclear magnetic resonance 
and photo-emission experiments. 

We now perform the minimization of the effective Hamil- 
tonian with NMTO set of parameters and compare the results 
with the equivalent minimization procedure performed in Ref. 
ifTlll with Mattheiss' parameters. As in Ref. 1,1 1.1 . we examine 
the following four magnetic phases: 

• AFM phase - all three in-plane bonds are antiferromag- 
netic; 

• RS phase - one in-plane bond is ferromagnetic and the 
other two are antiferromagnetic; (this is the spin struc- 
ture experimentally observed in V2O3); 

• ARS phase - one in-plane bond is antiferromagnetic 
and the other two are ferromagnetic; 

• FM phase - all three in-plane bonds are ferromagnetic. 

Fig[T]shows the plot of energy per V-ion, Ey, obtained with 
molecular (a) and atomic (b) variational functions using the 
NMTO set of parameters, as a function of Hund's coupling 
J and for U2 = 2.5eV. A direct comparison between Fig. [T] 
(a) and Fig. [T] (b) shows that for all values of J the energy 
of the atomic state is lower than the energy of the molecular 
state. This confirms what we had more qualitatively suggested 
above - that the molecular state is not supported by the NMTO 
calculations. 

The region of stability of experimentally observed RS mag- 
netic structure is again rather small, as we have obtained pre- 
viously with Mattheiss parameters for the molecular state. 
In this case, this region also shifts towards a lower value of 
Hund's coupling, J ~ 0.54 0.62eV. This value is close 
to J = 0.5 eV, extracted from optical studies of Qazilbash 
et al}^, as well as to the old value J — 0.59 eV given by 
Tanabe and Suganc^*^ by fitting optical spectra and at odds 
with the much bigger values ( J ^ 1 eV) recently used in the 
literature.'" 



For the orbital structure, the minimization gives the values 
of the orbital mixing angles 0^ ~ and ijji ~ 7r/2 at all 
sites for the RS configurations, ie, an orbital wave function 
iV'f) = |0)i = legCg). This is the solution without orbital 
degeneracy as found by Ezhov et al.}^ The absence of aig 
electrons in the solution can be understood as follows: the 
small gain in the kinetic energy, due to the significant reduc- 
tion of the hopping matrix elements in the NMTO calcula- 
tions, can not compensate the effect of the trigonal splitting, 
which pushes aig orbital to higher energies and favors the oc- 
cupancy of the |0) states on all the atoms. As a result, the 
ground state is a non-degenerate e^-doublet, and the orbital 
degeneracy is completely lifted, contrary to the ground states 
found in Ref. ll9lfTlll . where the important aig contribution to 
the kinetic energy in the molecule due to the aig — aig hop- 
ping term p allowed a finite occupancy of aig orbitals. It is 
worth mentioning that even with Mattheiss set of parameters 
it was possible to find a ground state le^Cg), however only for 
larger values of the trigonal splitting At > 0.4 eV. 

One more remark: due to the competing presence of the 
AFM, ARS and FM phases, the stabilization energy of the RS 
atomic solution is rather small, the first excited state is located 
about 2 meV above the ground state. Despite all differences 
between the ground state solutions obtained with NMTO and 
Mattheiss parameters, this is a common feature which is ap- 
parent from FiglUand from Figs. 5 and 8 of Ref. 1 1 1 ] : both the 
stability region and the stabilization energy for the the ground 
state with RS magnetic structure are much reduced in contrast 
with the spin 5=1/2 case considered by CNR^ . 

To illustrate more completely the properties of the effective 
spin-orbital model, Eq. |5] in Fig. |2]we present its ground-state 
phase diagram, in the parameter space specified by a/r and 
J /U2- All other parameters are taken from the NMTO set and 
we consider the trigonal splitting Af = 0.27 eV. ^From Fig|2] 
it can be seen that the ground state is successively changed 
from AFM to FM magnetic phase as the Hund's coupling 
increased. In order to gain maximum energy from orbital- 
dependent exchange terms, more complicated magnetic struc- 
ture are realized at the intermediate values of J, and that 
RS structure is stabilized in a strip between AFM and ARS 
phases. 

The examination of orbital structure for various ground 
states, presented on the phase diagram (Fig|2|i, has shown that 
orbital configurations in RS and ARS states are basically the 
same: the orbital degeneracy is lifted and all states are occu- 
pied by only Cg electrons. For the AFM and FM phases we 
find a continuum of orbital degeneracies with any mixing an- 
gles 9i and ■4>i, depending on value of the ratio J/C/2 and a/r. 

This phase diagram is very different from other phase dia- 
grams, which have been proposed for AFI of V2O3, i.e. from 
that of CNR model, ^ from the molecular model of Shiina et 
al}'' from that proposed by us in Ref. ifTTIl . The phase space 
occupied by RS phase is larger, and it happens at the realistic 
intermediate values of J jU^- 



J/U, 



AFM 

FM 




J 



FIG. 1 ; The energy per V atom as a function of J for different spin 
configurations (AFM, RS, ARS and FM type). The hopping parame- 
ters are defined by NMTO method'\ trigonal sphtting At = 0.27eV. 
In panel (a) the minimization is performed with molecular trial wave 
function while in panel (b) the minimization is performed with 
atomic trial wave function l[7]l. (Color online) 



IV. DISCUSSION 

In this paper we have studied the ground state of the effec- 
tive spin-orbital Hamiltonian (Eq. |2]l with the NMTO set of 
TB parameters. We have shown that the formation of stable 
molecules throughout the crystal is not favorable in this case. 
For all values of the Hund's coupling the minimization proce- 
dure with an atomic wave function gives a lower energy than 
in the minimization with a molecular wave function. The anal- 
ysis of the orbital structure of the ground state shows that the 
orbital degeneracy is lifted, and that two t2g electrons occupy 
the degenerate eg doublet. We also computed the mean-field 
phase diagram and showed that the RS structure is realized 
at realistic intermediate values of the Hund's coupling. The 
phase space of the RS phase is significantly enlarged com- 
pared to our previous study," although the gain of energy 



FIG. 2: Phase diagram in the (a/r, J/U2) parameter plane. Other 
hopping parameters are taken from NMTO set from Table I. Here 
AFM, FM, RS and ARS denote the corresponding type of magnetic 
order. Solid lines indicate the phase boundaries. (Color online) 



compared to competing phases is again small, of the order 
of ^ 2 meV. However, in every other aspect, the variational 
orbital ground state with two t2g electrons occupying the de- 
generate Eg doublet is in disagreement with the experimental 
data. 

First, it contradicts the observation of unquenched angu- 
lar momentum (< L >^ — O.S/is) found by Paolasini et 
al^"^ through the use of non-resonant magnetic x-ray scat- 
tering. Since < L >= over a real ground state, we 
need to introduce the spin-orbit coupling in the configuration 
jegEg >, obtaining the new ground state |e^e~ where 
= "^(TEg — *eg). Unfortunately this state still gives a 
zero expectation value for the orbital moment, due to its in- 
variance under complex conjugation (modulo an overall sign). 

Moreover it is now well establishe d'^'^^'^"^ that the resonant 
(1, 1, l)„i monoclinic reflection observed by Paolasini et al^^ 
is of magnetic origin. Specifically, it is due to a mixture of 
contributions of magnetic octupole and magnetic quadrupole 
moments of the V atom. The octupole moment results from a 
quadrupole-quadrupole {E2E2) x-ray transition of even par- 
ity, corresponding to the operator [[L (g) r] (X) f]^, whereas the 
quadrupole moment results from a dipole-quadrupole {E1E2) 
x-ray transition of odd parity, related to the spherical tensor 
[L r]^. However, the resonant scattering amplitude (which 
is related to the average values of both these operators in the 
|e^e~) ground state) vanishes, again due to the invariance of 
this state under complex conjugation. 

Yet another discrepancy comes from the fact that the orbital 
occupation of the \e'^eg) state disagrees with the population 
analysis of the t2g orbitals derived by Park et al^ This latter 
was based on polarized x-ray absorption measurements at the 
^2,3 edge of Vanadium in V2O3 and provided a value of 17% 
for the aig orbital occupancy in the AFI phase. 

Finally, a large value of the trigonal splitting A( together 
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with the same structure of ground state (le^eg >), as found 
in the paramagnetic phase'^ by the same NMTO implemen- 
tation of the LDA+DMFT approach, would lead to a highly 
anisotropic NMR line shape, with a sizable quadrupole split- 
ting in the V nucleus, and to a highly anisotropic Van Vleck 
susceptibility. 

In contrast, Rubistein^^ found, back in 1970, that in the 
paramagnetic phase of V2O3 there is no detectable anisotropy 
of the Knight shift and no observable quadrupole splitting of 
the NMR spectrum. Even earlier, Jones^^ performed suscep- 
tibility measurements for various crystalline orientations with 
respect to the magnetic field direction and found very little 
anisotropy. 

We emphasize instead that all the experimental observa- 
tions mentioned above are well explained in the framework 
of the molecular ground state picture^ ", which works espe- 
cially well after the inclusion of the spin-orbit interaction and 
the coupling with the lattice degrees of freedom'^. 

We therefore strongly doubt the applicability of the NMTO 
method to describe the physics of V2O3, even though the 
method in itself is quite sophisticated for band structure 
calculations and is certainly in keeping with the localized, 
Hubbard-type of approach to the description of correlated sys- 
tems. 

In fact we note that the difficulties found in the framework 
of the NMTO-method when applied to V2O3 reflect a gen- 
eral problem of all LDAh-U and LDAh-DMFT calculations 
based on local atomic correlations inside a single Vanadium 
atom. All these approaches lead essentially to the same kind 
of ground state, with two eg electrons in spin one configura- 
tion without orbital ordering and an almost empty aig band. A 
similar result was found by CNR in their second paper of the 
series^, starting with an empirical band structure and an unre- 
stricted Hartree-Fock calculation, which is very similar to the 



modern LDAh-U scheme. 

In a previous publication^^, one of us has already high- 
lighted the main drawbacks of implementation of DFT-based 
ab-initio methods to V2O3. This compound is a very peculiar 
strongly correlated system, in which non-local correlations 
play an important role. The molecular ground state of V2O3 
is stabilized because of these non-local interactions (e.g., cor- 
related hopping along vertical bond), the physics of which can 
not be captured within a method based on local interactions. 
The NMTO-downfolding belongs to the same class of DFT ab 
initio methods and, therefore, holds all shortcomings of local 
correlated approaches. Stated differently, non-local correla- 
tions would lead to non-local one particle effective potentials 
that might substantially modify the shape and the range of the 
Wannier-like functions. 

One step in this direction has been provided by a new gener- 
ation of correlated non-local quantum chemical ab initio cal- 
culations. In this approach a combined exact diagonalization- 
ab initio method (EDABI) by Spalek et afi^ is used as a 
promising method to estimate correlated hopping. This latter 
is based on a definite procedure of construction of the many- 
particle trial wave function expressed in terms of an adjustable 
one particle basis, followed by the solution of a self-adjusted 
non-local and non-linear wave equation obeyed by the basis 
functions. We believe that the EDABI method is particularly 
well suited for treating situations in which the inter-particle 
correlations are not weak and lead to non-local interactions. 

Still, the major outcome of the NMTO-downfolding 
method is the realization that in-plane hopping parameters are 
important and a model of weakly interacting vertical corre- 
lated molecules, which neglects in-plane hopping, is grossly 
oversimplified. The interactions between molecules, originat- 
ing from in-plane hopping, are therefore not negligible and 
should be included into the theory. 
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